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Abstract: In the presence of dense matter quantum anomalies give rise to two 
new transport phenomena. An anomalous current is generated either by an external 
magnetic field or through vortices in the fiuid carrying the anomalous charge. The 
associated transport coefficients are the anomalous magnetic and vortical conduc- 
tivities. Whereas a Kubo formula for the anomalous magnetic conductivity is well 
known we develop a new Kubo type formula that allows the calculation of the vor- 
tical conductivity through a two point function of the anomalous current and the 
energy current. We also point out that the anomalous vortical conductivity can be 
understood as a response to a gravitomagnetic field. We apply these Kubo formulas 
to a simple holographic system, the R-charged black hole. 
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1. Introduction 

Anomaly induced charge separation effects in heavy ion collisions are currently under 
intense discussion [1]. There are basically two such effects called the chiral magnetic 
effect and the chiral vortical effect. Both effects rely essentially on the presence of 
an imbalance between left-handed and right-handed quarks. In many models this 
imbalance is described by introducing a chemical potential for the axial charge. The 
chiral magnetic effect describes the generation of a current in a background magnetic 
field and the chiral vortical effect describes the generation of a current in the presence 
of a vorticity field in the charged fluid. 

The chiral magnetic effect has actually a somewhat complicated history as it 
seemingly has been discovered independently several times [2,3]. In any case its rel- 
evance to heavy ion collisions has been argued for in [4] and from then on there has 
been an intensified activity in deriving and understanding it from various perspec- 
tives. Most of the work has been done from a traditional field theory approach [5-18]. 
Since it is related to the anomaly one would naturally believe that the chiral mag- 
netic conductivity is an exact result and should apply also at strong coupling. Indeed 



- 1 - 



it has recently been shown that higher loop corrections do not contribute and thus 
the result is indeed exact in perturbation theory [19]. Non- anomalous corrections to 
it due to interactions have been however argued for in [20]. Lattice calculations of 
the chiral magnetic effect have been reported in [21-23]. Recently it also has been 
shown that at finite magnetic field there are now propagating hydrodynamic modes 
called chiral magnetic waves [24] . 

Results reported in holographic models were however consistent with a vanishing 
chiral magnetic conductivity or used the consistent form of the anomaly with a non- 
conserved vector current [25-27] ^. This discrepancy in the holographic calculations 
has been resolved in [29]. There it was pointed out that a chemical potential can 
be introduced in various gauge equivalent forms. Since gauge invariance is however 
violated by the anomaly these different formalism have to be kept distinguishable 
during the calculation and it was shown that the weak coupling result corresponds 
to a formulation in which the chemical potential is not represented by a background 
value for the temporal component of a fiducial gauge field but rather by certain 
boundary conditions on the fields. In fact this seems the most natural way of intro- 
ducing a chemical potential for an anomalous gauge symmetry, since after all vector 
fields that couple to anomalous currents are absent in nature.^ 

The anomalous vortical effect on the other hand has first been noticed in a holo- 
graphic calculation of the hydrodynamics of R-charged asymptotically anti de Sitter 
black holes [31,32]. Soon it was pointed out that such a term in the constitutive 
relations for an anomalous current is not only allowed but even necessary in order 
to obtain an entropy current whose divergence is positive definite [33]. Recently it 
was suggested that the interplay of chiral magnetic effect and chiral vortical effect 
leads to interesting signatures of baryon number separation and charge separation in 
heavy ion collisions [34] whereas [35] argued for enhanced production of spin-excited 
hadrons. The conserved electromagnetic current has been explicitly incorporated in 
the hydrodynamic constitutive relations recently on [36] whereas [37, 38] considered 
an extension with global or gauged abelian and non-abelian symmetries. Very re- 
cently a STU model realizing the chiral magnetic and the chiral vortical effects has 
been presented in [39]. 

Let us briefly discuss how a chemical potential can be introduced in fleld theory. 
Often it is viewed as a deformation of the theory by a new coupling of the form 

^See also [28]. 

^In [30] it was argued that a chemical potential should be introduced only for non-anomalous 
charges. By addition of a Chern-Simons current a truly conserved axial charge can indeed be 
defined. This charge is however only gauge invariant if it is integrated over all of space. Therefore 
this approach seems not directly applicable to situations where the charge is supposed to be localized 
inside the small fireball of heavy ion collision. We rather prefer to work with chemical potentials 
for anomalous charges that are not truly conserved but localizable within small regions in a gauge 
invariant way. 
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H ^ H — jiQ where H denotes the Hamiltonian, jj, the chemical potential and Q 
the charge operator. If Q is the zero component of a conserved current then we 
can think of this as arising from coupling the current to a (fiducial) gauge field 
and choosing a constant background value for the temporal component of the gauge 
field Aq — ji. In a finite temperature context fields obey the KMS state condition 
(t){t + i/T) — r](j){t), where 77 = ±1 for either Bosons or Fermions. Due to the 
underlying gauge invariance we can however also introduce the chemical potential 
in a different but equivalent way. We can perform a gauge transformation that sets 
Ao = 0. That does of course not mean that the chemical potential has vanished from 
our theory. Since the (imaginary) time direction is compactified due to the KMS 
condition such a large gauge transformation changes the boundary conditions on the 
fields. In fact the correct KMS state conditions in the gauge with A^ = are 



where q is the charge of the field (p. It seems useful to think of the first formulation as 
a deformation of the dynamics of the system whereas the second formulation leaves 
the dynamics, i.e. the Hamiltonian H unchanged but modifies the state space (see 



If we try now to define an analogue of the chemical potential in the case of 
an anomalous theory we have to be aware of the fact that these two formulations 
might give different results since a gauge transformation on the fields is not anymore 
an invariance of the quantum theory. Following [29] our approach in holography is 
to choose a formulation that distinguishes between the elementary definition of the 
chemical potential as the energy that is needed to introduce a unit of charge into the 
system from the actual value of the temporal gauge field. In the holographic setup the 
energy needed to introduce a unit charge into the system is given by the potential 
difference between the horizon and the boundary. This quantity does not change 
under gauge transformations even in the anomalous case. This condition describes 
the state of the system. On the other hand we can then allow for an arbitrary gauge 
field background that is not a priori related to the chemical potential. Choosing 
such a non- vanishing value for will then also deform the dynamics. In [29] it was 
shown that calculating the chiral magnetic conductivity with the deformed dynamics 
leads to a zero result but that a non zero result coinciding with the weak coupling 
result [6] is obtained if this deformation is switch off. 

Whereas the actual chiral magnetic effect relevant to the physics of heavy ion 
collision depends on the interplay of a conserved vector current and a non-conserved 
chiral U{1) current we will focus in this paper on the most simple setup where only 
one anomalous U{1) current is present. One of the main results of this paper is a 
Kubo formula for the anomalous magnetic and vortical conductivities ^b, iv as they 
are defined through the constitutive relations in [31-33]. These Kubo type formulas 
are: 




e.g. [40-42]). 
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Here are the spatial components of the anomalous current, T*^ is the energy flux 
and kn the momentum. The correlators are retarded ones. 

That we have to compute the sum of two correlators has to do with the particular 
way the coefficients are defined in the constitutive relation and is related to an 
ambiguity in the definition of the fiuid velocity at first order in derivatives. We also 
have explicitly included in our definition of the anomalous magnetic conductivity 
that the temporal component of the gauge has to vanish. It is remarkable that the 
anomalous vortical conductivity does not depend on the gauge field background and 
therefore can be evaluated for any Aq. 

This paper is organized as follows. In section two we derive the Kubo formulas 
for the anomalous conductivities from the constitutive relations derived in [33]. We 
point to an inherent ambiguity in the definition of the fiuid velocity. Furthermore we 
define new elementary conductivities relying on the interpretation of the anomalous 
vortical effect as a gravitomagnetic effect and relate them with the Landau frame 
conductivities defined above. 

We then go on and define our holographic model, which is nothing but (a trunca- 
tion of) the asymptotically R-charged black hole background. We define the current 
and show that it obeys an anomalous conservation law. We compute then the rele- 
vant matrix of retarded Green functions analytically in the hydrodynamic limit and 
show that applying the Kubo formulas (1.2), (1.3) the results established in [33] are 
obtained. We find that the constitutive relations for the anomalous current have to 
be modified by including a Chern-Simons coupling. We would also like to comment 
that the computation of the retarded correlators for the holographic model consid- 
ered in this work has been addressed before in [43], though the ambiguity in the 
fiuid velocity was not addressed and only the magnetic conductivity was computed. 
In fact, the Green functions found by the authors correspond to the particular case 
of setting that velocity to zero, i. e. neglecting the energy fiow due to the external 
background magnetic and vorticity fields^. We then proceed and compute frequency 
dependent conductivities numerically. 

In section four we draw our conclusions and discuss possible further applications 
of the Kubo formalism. In particular we emphasize that the Kubo formula allows 

^Let us also note that chiral dispersion relations in the shear sector have been found in [44]. 
This is however a higher order effect and not captured by first order hydrodyamics which is our 
prime interest here. 
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for a first principle calculation of the chiral vortical conductivity along the lines of 
the calculation of the chiral magnetic conductivity of [6,45]. 

In Appendix A and B we collect the lengthy expressions for the perturbed ac- 
tion and the analytic solutions for the fluctuations in the hydrodynamic regime. In 
Appendix C we explore the dependence of the correlators and the conductivities in 
the fluid velocity at zero frequency. 

2. Kubo formulas and constitutive relations 

Hydrodynamics can be understood as an effective, dissipative field theory for a fluid. 
As is usual for effective field theories it is organized as a derivative expansion. To 
lowest order one starts with the equilibrium expressions of energy momentum tensor 
and currents. Dissipative effects appear at first order in derivatives. According to the 
underlying space time symmetries the derivative terms can be organized in certain 
tensor structures. The most important input parameters to hydrodynamics are then 
the coefficients in front of these derivative tensor structures. These are the transport 
coefficients, such as viscosities or conductivities. These constitutive relations can 
also be written down in the presence of generic external fields, such as electric and 
magnetic fields or gravitational fields. 

The constitutive relations for a relativistic fluid carrying charge of an anomalous 
f/(l) symmetry have first been obtained via holographic methods in [31,32]. They 
have been further studied from a purely hydrodynamic point of view in [33]. These 
constitutive relations are 



Here e is the energy density, P the pressure, n the charge density and the local 
fluid velocity. The fluid velocity can be chosen such that the dissipative parts t^'^ 
and u'^ obey u^t^'^ = and u^u'^ = 0. 

In the presence of an external background gauge field and assuming the current 
to be anomalous the conservation equations take the form 
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where P'"' = gf"" + u/'u", = F'^'u^, 5^ = le^^^^^M^F^A and uf" = e^^'P^u^dpUx is 
the vorticity of the fluid. In addition to the usual transport coefficients, the shear 
viscosity 77, the bulk viscosity C and the conductivity a there are two new, parity 
violating terms proportional to the magnetic field and the fluid vorticity These 
are induced by the anomaly. 

We now want to derive Kubo formulas for the anomalous transport coefficients. 
First let us assume that the fluid is close to rest = {l,v^,Q,v^) with only small 
velocities such that u'^u^ = — 1 + 0(2). If we interpret the constitutive relations 
as the one point functions of the current and the energy momentum tensor in the 
presence of external sources, we can obtain the two point function of currents by 
differentiating with respect to the gauge potential and the metric^. 

In the constitutive relations there are no terms that depend on the metric only. 
This is a consequence of the fact that covariant curvature expressions are second 
order in derivatives and therefore can not appear as long as we restrict ourselves to 
first order hydrodynamics. Coupling to the metric is achieved simply by taking all 
derivatives as covariant ones. 

We now assume that in addition to the velocity field we have a metric pertur- 
bation of the form g^j^i, — -q^i, + h^j^^ with only htx, hfz, hy^ and hy^ different from 
zero. The gauge field perturbations we take to be — (0, a^;, 0, a^). All unknown 
functions are supposed to depend only on {t,y). To first order in derivatives, sources 
and the velocity the relevant constitutive relations are now 

jm ^ ^^m _ ^mn^^ ^g^j^^^ ^ g^^n^ _ ^g^^^ _ ^^^^nng^^^ ^ ^^2) , (2.7) 
rj^m _ ^ p^^m ^ pj^^^^ ^ 0(2) , (2.8) 

= -Phy^ - l^{dyV^ + dthyra) + 0(2) , (2.9) 

with m e {x, z}. In addition we have the (non-) conservation law 

[(e + P)dt - vd^y]v"' + ndtam + (e + P)dthtm - vdtdyhym = 0(3) . (2.10) 

The orders on the right indicate that these equations hold up to n-derivatives. We 
could now eliminate the velocity from the conservation equation of the stress tensor. 
However, this equation does have a non-vanishing zero mode to first order in the 
derivatives. For simplicity let us take the zero frequency sector. We see then that any 
velocity of order 0{dy) solves the conservation equation. Even at non-zero frequency 
we could still have a velocity of the form v.^ ~ (e+p)af-r?92 /(^i v) fo^^ an arbitrary 
function f{t,y). In non-anomalous theories we can exclude such terms since parity 

^Including a fiducial metric in the constitutive relations of hydrodynamics for the purpose of 
deriving Green-Kubo-Mori formulas goes back to Luttinger's theory thermal transport coefficients 
[46], where he states that if the gravitational field didn't exist, one could invent it for the purpose 
of deriving the transport coefficients. 
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conservation together with the 50(3) rotational symmetry would tell us that the 
spatial momentum can enter only as k'^ in the equations. In the presence of the parity 
violating anomaly we can not assume this anymore and have to take the kernel of 
the energy-momentum conservation into account. We will denote this kernel as Vm^ 
and treat it as an arbitrary function of the sources. 

Let us further simplify to zero frequency. We find then 



rptm 

jm 



(e + P)t;(i) + Ph 



tm 5 



nv. 



(1) 



(2.11) 
(2.12) 



We can differentiate now with respect to htm or an and eliminate the undetermined 



VmT^i and (vmJ^) correlators and arrive at 
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5'"" + eye™"'ifci + 0(F).(2.14) 



From this the Kubo formulas (1-3) and (1.2) follow immediately. 

In fact even without using the (non-) conservation equation for the energy mo- 
mentum tensor we can find some useful relationships between Green functions. We 
can solve the constitutive relation (2.8) for v"^ and plug it into the one for the cur- 
rent. Upon differentiating with respect to the metric component hm or the gauge 
field an we find that for m = x, n = 2; the following relations hold 



P) 



ky e + {T^^T*') 



rr) 

JXrptz^ _ 



n 

e + P 
n 



p 



{T'r) + 0{k), (2.15) 



P 



^rptXrptZ^ 



--^v + 0{k). 
P and (r*V^~) 



(2.16) 



and we recover 



To lowest order in uj and k we have ^T*^T*^) 
the previous Kubo type formulas. 

Finally we need to understand why the zero- mode velocity Vm^ is present. We 
know that a variation in the charge distribution will cost us an energy of the form 
6e = fi6Q. If we imagine a test charge 6Q moving through the charged plasma 
it will therefore generate a current 6 J and induce also an energy current of the 
form = fi6J^. For a finite current we should integrate this over /i and obtain 
T** = fi'^d^' . In the Landau frame this implies however the generation of a 
velocity through T*' = (e + P)w*. If we are interested in the charge current that sits 
on top of the one produced by the fluid velocity we should subtract the component 
of the current that is due to the drag of the fluid in the Landau frame 

n 



P 

drag 



nv 



e + P' 



(2.17) 
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From what we outhned above it follows then that the transport coefficients 
and that appear in the Landau frame constitutive relations are^ 

= = <Jb —B / /i'd(7B(/i') , (2.18) 

Cv^^av- -— / li'davili') . (2.19) 

Using these relations we find the following form of the relevant correlators 

(J™ J") = ie^^^^ki ae , (2.20) 

(r*"*J") = ie^'^iki [ n'daeifJ.') , (2.21) 
Jo 

(T*™T*"> = ie™"'A;z / n'^dasilJ.') ■ (2.22) 
Jo 

The second line in these equations fixes also the chiral vortical conductivity by 
complex conjugation 

^jmj:tn^ ^ ^ytnjm^t = ie'^^.l ki(Tv ■ (2.23) 

We also would like to point out now that the vortical effect can be understood 
as a anomalous gravitomagnetic effect. For slowly rotating and well localized matter 
distributions it is often convenient to write the gravitational field as 

ds^ = -(1 - 2^g)dt^ - 2Agdxdt + (1 + 2^g)5ijdx'dx^ . (2.24) 

$g is the Newtonian gravitational potential whereas Ag is the gravitomagnetic po- 
tential. In full analogy we can define the new elementary chiral magnetic and vortical 
conductivity through 

J^aeB, (2.25) 
J^avBg, (2.26) 

where we introduced a "gravitomagnetic" field Bg = e^^^djAg^^ for the metric fluctu- 
ation®. 

We will now solve for the retarded two-point functions in the holographic model. 
We do this analytically in the hydrodynamic limit. Here we will find an analogous 



^We would also like to point to [47] where using an effective field theory approach it has been 
argued that the terms of higher order in the chemical potential of the anomalous conductivities are 
ambiguous due to infrared effects. 

^Scc [48] for a review on the gravitoelectromagnetic approximation to general relativity. 
The gravitoelectric field Eg = — Vi>g has of course also a nontrivial effect and measures the 
(thermo)electric conductivity which is most easily seen when we substitute the temperature in 
the constitutive relations by the gravitational potential dependent local Tolman temperature [46]. 
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ambiguity that appears as an undetermined integration constant and that is related 
with the arbitrariness in the definition of the fiuid velocity already mentioned. This 
ambiguity disappears once the correct boundary conditions are imposed. Numerically 
we can work to all orders in frequency and momentum and there we will find that 
all correlators are uniquely determined. 



3. Kubo formulas in Holography 

The holographic dual description of a 4-dimensional U{1) chiral gauge theory is 
simply given by a 5-dimensional Einstein-Maxwell model supplemented with a topo- 
logical term ^, 

S-yIgI "^'"^^ + " i^' + ^e'^^'^^AMF^nFcn^ , (3.1) 

where we can keep the value of the Chern-Simons parameter k to be arbitrary. Latin 
indices denote bulk coordinates whereas Greek indices denote boundary coordinates. 
The equations of motion 



= Rab-^{R+ji^\9ab + ^FacF^ b + ^F'qab , 



= V^F^^ + «:e^^^^^FMivFcD , (3.2) 
admit the following exact AdS Reissner-Nordstrom black-brane solution 

A = (l){r)dt={/3-^)dt, (3.3) 

where the horizon of the black hole is located at r = th and the blackening factor of 
the metric is 

J.4: ^6 

The parameters M and Q of the RN black hole are related to the chemical potential 
II and the horizon rn by 



/('■) = i--;5- + ^- (3-4) 



The Hawking temperature is given in terms of these black hole parameters as 



^For the sake of clarity let us mention that here we are defining the e tensor to be related with 
the Levi-Civita symbol by eABCDE = y^e{ABCDE), with e(0rf23) = e(0123) = -1. This 
normalization is consistent with having positive Chern-Simons parameter k. 
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The pressure of the dual gauge theory is 



The theory is conformal obeying e = 3P. 

Accordingly to the discussion in the introduction, we identify /i as the chemical 
potential of the gauge theory since it is precisely the gauge invariant quantity corre- 
sponding to the energy needed to introduce a unit charge into the system, i. e. it is 
the difference between the scalar potential at the boundary and at the horizon. On 
the other hand, the parameter /3 corresponds to the boundary value of Aq and as 
already mentioned, it can be viewed as a deformation of the dynamics of the system. 
Following the arguments in [29], it will be eventually set to zero in order to recover 
the weak coupling result for the chiral magnetic conductivity, but for the moment 
we keep it as an arbitrary gauge background for the boundary theory. 

The study of the effect of the anomaly in the hydrodynamic behavior of the 
chiral gauge theory and the extraction of the corresponding transport coefficients 
passes through computing the holographic two-point correlation functions of the 
current and the stress tensor. For doing so, we expand the action to second order in 
perturbations of both the gauge field and the metric. 

Am Am + aM ] Qab Qab + hAB ■ (3.8) 

As a first step, we use the action to first order in perturbations to define the 
current coupled to the gauge field and the stress-energy tensor of the dual field theory. 
To first order in gauge fluctuations, the action on-shell reduces to a boundary term 

-TtalH-^' ^ f ^--"^.fc.) \_ . (3.9, 

from which we can read off the boundary current, that applying the holographic 
dictionary is given by 



6S fr^ur 4fi: 



dA^{r — )■ oo) IGnG \ 3 



(3.10) 



The charge density can be computed as the time component of the current and is 

It can be easily checked that this current is not conserved. Using the equation of 
motion for the gauge field in (3.2), the divergence of the anomalous current reads 



E^B^, (3.12) 



67rG 
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where E^^ and correspond to the electric and magnetic fields of the boundary 
theory. The result above differs from that of [33] in a numerical factor that can be 
explained taking into account the different normalization of the gauge field on top 
of the fact that we are using different definitions of the current: whereas we define it 
through its source term in the action, the current in [33] is defined as the subleading 
term in the near boundary expansion of the gauge field ^. Taking these differences 
into account we find that our anomaly coefficient c is 1/3 of the anomaly in [33]. 

In an analogous way, we can work out the stress-energy tensor of the gauge 
theory from the bulk action to first order in perturbations of the metric. In this 
case, it is necessary to include the Gibbons-Hawking term in order to keep a well- 
defined variational principle and also a counterterm coming from regularization of 
the boundary action [49], 

where 7 is the induced metric on the AdS boundary, K is the extrinsic curvature and 
i?(4) is the Ricci scalar of the boundary metric. The Chern-Simons term does not 
depend on the metric since it is a topological contribution, hence the stress tensor 
will not receive any correction due to its presence. In fact, one recovers the result 
for the energy-momentum tensor obtained in [49] from a generic five dimensional 
asymptotically AdS geometry. 

The action to second order in perturbations of the metric and the gauge field 
will receive contributions from the bulk action (3.1), from the Gibbons-Hawking 
term and from the gravity counterterm (3.13), but also from the counterterm coming 
from regularization of the gauge part of the action needed to cancel the logarithmic 
divergence. This term is given by 

= -TTTT^ log^ / d'x^F.^F'^-^ . (3.14) 
IottG J 

The precise form of the contribution of each of these terms to the second order action 
can be found in Appendix A. 

Notice that the expansion of the action includes both terms to second order in 
gauge fluctuations and in metric fluctuations, but also mixed terms. In principle one 
could think that only the mixed terms will contribute to the off-diagonal correlators, 
but of course this is not true due to the holographic operator mixing under the renor- 
malization group flow. Thus a priori all the terms contribute to all the correlation 

8 The precise relation is J/,;^.^ = i (-^sonSurowka + 3iM^''''''''^'^^^P>) ^^^""^ Sauge fields 
and Chern-Simons coupling differ by a factor of | from [33]. 
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functions. Collecting the perturbation fields om and Hmn in a single vector $/ and 
inserting its Fourier mode decomposition, 



^\r,xn = J 0^,<{r)e-''-'^''' , (3.15) 

the complete second order action on-shell can be compactly written as a boundary 
term 

^^^'^ = I Tfr^i'^-^^^^'^'/ + ^^-k^iM) , (3.16) 

J [Its:)'' r^oo 

where derivatives are taken with respect to the radial coordinate. In order to avoid 
double counting and to keep the correct contact terms in the Green functions one 
should be careful in the extraction of the contribution to the B matrix coming from 
the bulk action. 

We can compute the holographic response functions from (3.16) by applying the 
prescription of [50,51]. As first noticed in [52], for a coupled system the holographic 
computation of the correlators amounts to find a maximal set of linearly independent 
solutions that satisfy infalling boundary conditions on the horizon and that source 
a single operator at the AdS boundary. Following [53], we can construct a matrix of 
solutions j{p, r) such that each of its columns corresponds to one of the indepen- 
dent solutions and that at the boundary becomes the unit matrix. Therefore, given 
a set of boundary values for the perturbations, (fil, the bulk solutions are 

^iir)=F^jik,r)^i. (3.17) 

The F matrix actually corresponds to the bulk-to-boundary propagator matrix for 
a coupled system of bulk fields, for which the dual field theory operators are well- 
defined at the UV scale and get mixed under RG flow. Inserting these solutions in 
the second order boundary action, the holographic Green functions are finally given 

by 

Guik) = -2 hm {AiMiF''j{k,r)y + Bij) . (3.18) 

Let us now present the precise setup that we will consider from now on in order 
to study the anomaly effects. Without loss of generality we consider perturbations 
of frequency u and momentum k in the y-direction. We discuss only the shear chan- 
nel (transverse momentum propagation) since the off-diagonal correlators appearing 
in the Kubo formulas for the chiral vortical and the chiral magnetic transport co- 
efficients belong to the vector channel. This implies that we have to switch on the 
fluctuations Ai, h\ and h^y, where i = x,z. For convenience we deflne new parameters 



(3.19) 



- 12 - 



and the compact coordinate u = r^/r"^, for which the horizon sits at u = 1 and the 
AdS boundary at m = 0. 

Finally we can write the system of differential equations for the shear sector, 
that consists on six second order equations and two constraints 

= coh'i{u) + kf{u)h';l:{u)-3auuBi{u), (3.20) 

^ .m.>^ O^'ftiM + ^^feiW) .3auBHu) , (3.21) 

u uf{u) 

= K\u)+(Q^--) Kiu)-^^ikujhliu)+u'Kiu)), (3.22) 



fiu) uJ ^ ufiu] 



/(m) uf{uY f{u} f{u) 

where we define Bi{u) = Ai{u)/n and K = AnL^fi/rf^. There still remains a residual 
gauge symmetry on the gravity sector 6hl = uX and Sh]. = —kX. This symmetry 
tells us that the metric components are not independent. In fact, it is possible to 
combine (3.20) with (3.21) to obtain (3.22). 

Since we are interested in computing correlators at zero frequency, we can drop 
out the frequency dependent parts in the equations and solve the system up to first 
order in k. In this limit, the fields h], decouple from the system and take a constant 
value. The reduced system can be written as 

= h'!\u) - - 3 a n B'(u) , (3.24) 

u 

- B, [u) + ^^^^ B^[u) I e., K ^^^^ ^^^^ , (3.25) 

where we expand the fields to first order in the dimensionless momentum p = k/A-nT 

hliu) = hf\u)+ph^\u), (3.26) 
B,iu) = Bf\u)+pBf\u). (3.27) 

In order for the solutions to be physically sensible, they have to satisfy certain 
boundary conditions. The first condition is that they source the desired single op- 
erators in the UV, so the bulk fields must satisfy h\{^) = hi, -Bj(O) = Bi, where 
the 'tilde' parameters are the sources of the boundary operators. The second condi- 
tion comes from imposing infalling boundary conditions at the horizon. But what is 
the 'infalling' condition at zero frequency? At first sight, the condition leads to an 
ambiguity, since both solutions for the metric fluctuations are perfectly regular, al- 
lowing us to take an arbitrary horizon value. However, if we analyze the near horizon 
behavior of the fields for arbitrary frequency, 

hi{u) ~ (1 - n)^^'"/^"^+i , (3.28) 
Si(M) ~ (1 -m)-'"/^"^, (3.29) 
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and then take the hmit a; — )■ 0, we see that infalhng in this case means: regular 
for the gauge field and vanishing for the metric fiuctuation. In the appendix B we 
present the solutions to the equations of motion at = that satisfy this boundary- 
problem. The apparent freedom in fixing that integration constant is related with 
the zero mode found for the fiuid velocity in first order hydrodynamics. We devote 
Appendix C to explore this relation and show that the transport coefficients C,b and 
are independent of this parameter as one would expect. 

Now that we have the solutions for the perturbations, we can go back to the 
formula (3.18) and compute the corresponding holographic Green functions. If we 
consider the vector of fields to be 

$^(n) = (S^n), /i^(n), B,{u), h]{u)^ , (3.30) 
the A and B matrices for that setup take the following form 

^ = T^^Ti Di^g «/' -3 -) > (3-31) 
IottGL^ \ u u ) 



B 



AdS+d 













—3a 

__3_ 















—3a 

/ 



(3.32) 



/O \ 











B 



CT 











° ^ J 



(3.33) 



where we have split the B = BAds+d + Bct matrix into the contribution coming from 
the bulk and the Gibbons-Hawking actions, BAds+d^ and the contribution coming 
from the counterterms, Bct- With these matrices and the perturbative solutions 
we can construct the matrix of propagators. The non-vanishing retarded correlation 
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functions at zero frequency are then 



Gx.te - G,^t, - ^^^3 , (3.34) 
__ iV3kQK ik(3K 

Gx.tZ = GtX,Z = ~Gz±X = —Gtz.X = A ^ A ) (3.36) 

M 

Gtx,tx = Gtz,tz = 77 — TTT^ 5 (3.37) 



167rGL3 



Gtx,tz — ~Gtz,tx — — —^ — T^m — • (3.38) 



H 

The off-diagonal current-current correlator (3.35) is consistent with the axial-axial 
current correlator computed on [29]. As argued in the introduction, a background 
gauge field corresponds to a deformation of the dynamics of the boundary theory. 
For an anomalous theory, this deformation would lead to the appearance of an extra 
Chern-Simons term in the current constitutive relation 

which is not gauge invariant but is allowed because of the anomaly. If we deform 
the theory turning on a constant background Aq = /3, we get a contribution to the 
magnetic conductivity coming from (3.39) which is precisely the /3-part on (3.35). 
This extra contribution to the constitutive relation is just the difference between 
the current defined as the subleading term in the boundary expansion of the gauge 
field and the definition we are using, i.e. as the variation of the action including the 
terms stemming from the Chern-Simons part. Therefore, a non-vanishing f3 can be 
interpreted as a source for such a topological coupling. In fact, differentiating with 
respect to this background field it is possible to compute the three point function of 
currents 

{J\k)J^{-k)J\0)) = -^eij = -zkcu,. (3.40) 

On the other hand, if we switch off the deformation, i.e. vanishing background 
field, we recover the conductivities (2.20), (2.23), (1.2) and (1.3) 



a. = ^ = (3c)/., (3.41) 



-v = ^ = (3c)^, (3.42) 



47rG4 ' ' 2 
Vsg (ML2 + 3 4) 



3Q'k 3c f , 2 n^^^ \ 
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The last two expressions agree precisely with [33] once we take into account the 
difference in the anomaly factor c and the fact that we also dropped a factor 1/2 in 
our definition of the vorticity. They also obey of course the relations (2.21), (2.22) 
and (2.23). The result for as coincides with the previous holographic computations 
of [26, 29] for the chiral magnetic effect in an axial anomalous theory, and therefore 
with the weak coupling field theoretical result of [6], whereas the chiral vortical 
conductivity ay is a new result. 

3.1 Frequency dependence 

In order to study the frequency dependence of the chiral conductivities, we can use 
(2.20) and (2.23) to define 

aBiuj) = hm -^e™, ( JV^) , (3.45) 
avioj) = lim ^e^i, ( JT*^) . (3.46) 

It is important to notice that these, and not the C,v and ^b, are the relevant 
conductivities at finite frequency. The latter correspond to the conductivities mea- 
sured in the local rest frame of the fluid, where one subtracts the contribution to 
the current due to the energy flux generated when we put the system in a back- 
ground magnetic or vorticity field. But as we have seen, there is an ambiguity in 
the definition of the local rest frame: the fluid velocity is frequency and momen- 
tum dependent, one can just define it order by order in the hydrodynamic gradient 
expansion up to an arbitrary contribution. This automatically implies that the 
and are only meaningful in the zero frequency limit. On the other hand, the ay 
and cb conductivities are not subject to this problem: they capture the complete 
response of the system to the external magnetic fields. Therefore, it is sensible to 
define the frequency dependent chiral conductivities above, in an analogous way as 
for the A.C. electric conductivity. 

To study that dependence holographically, we have to resort to numerics. The 
nature of the system allows us to integrate from the horizon out to the boundary, 
so we should fix boundary conditions at the first one, even though we would like 
to be free to fix the AdS boundary values of the fields, hence the operator sources. 
Imposing infalling boundary conditions, the fluctuations can be written as 

hliu) = {1 - u^-^' Hliu) , (3.47) 
Kiu) = {1 - Hi,{u) , (3.48) 
B\u) = {l-u)-''"b\u) , (3.49) 

where w = uj/4ttT. As we saw, the remaining gauge symmetry acting on the shear 
channel implies that h] and h], are not independent. So if we fix the horizon value 
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of the {6*, HI} fields, the constraints (3.20) fixes 

3m6'(l) + {i + w)Hi{l] 



(2 - a)p 



(3.50) 



In order to find a maximal set of linearly independent solutions, we can construct four 
of them using linearly independent combinations of these horizon free parameters. 
In this way we construct the following independent horizon value vectors 



1 \ 




/ 







1 


3ia 






{2-a)p 




{2-a)p 



















V / 



V / 



V 







1 



3ia 



(2-a)p 



I 



\ 








1 

i+w 



\ 



(3.51) 



{2-a)p 



J 



The remaining two are given by pure gauge solutions arising from gauge transforma- 
tions of the trivial one. We choose them to be 



$(m) 



{ ° ^ 


/ \ 


w 





-p 














w 


^ / 


\-p) 



(3.52) 



Using the corresponding solutions we construct the F matrix of (3.17) in this way: 



(3.53) 



where Hj{u 



In Figure 1 is illustrated the behavior of the vortical and magnetic conductivities 
as a function of frequency for two very different values of the dimensionless temper- 
ature r = 27Tr}iT/fi. Both of them go to their corresponding zero frequency analytic 
result in the w — ?■ limit. The frequency dependent chiral magnetic conductivity was 
also computed in [26], though in that case the possible contributions coming from 
metric fluctuations were neglected. Our result for (Tb(w) agrees pretty well with the 
result found in that work in the case of high temperature when the metrics fluctu- 
ations can be neglected, but it develops a dip close to w = when temperature is 
decreased (see Figure 2), due to the energy flow effect. For small temperatures, the 
chiral magnetic conductivity drops to ~ 1/3 of its zero frequency value as soon as 
we move to finite frequency. The behavior of ay is slightly different: the damping is 
much faster and the imaginary part remains small compared with the zero frequency 
value. 
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Figure 1: Chiral vortical (up) and magnetic (bottom) conductivities as function of the 
frequency at r = 36.5 (left) and r = 0.24 (right). Red doted points represent real part and 
thick blue line the imaginary conductivity. 




0.005 0.010 0.015 0.020 4;rr 0.005 0.010 0.015 0.020 4 irT 



Figure 2: Chiral vortical (up) and magnetic (bottom) conductivities as function of the 
frequency close to w = 0. Real (left) and imaginary (right) part of the normalized conduc- 
tivity for different values of the dimensionless temperature. 



In Figure 2 we made a zoom to smaller frequencies in order to see the structure of 
the dip on aB and the faster damping on ay- In Figure 3 we show the conductivities 
for very small temperature. From this plots we can infer that at zero temperature 
the conductivities behave like as = aa^ (l + ^-^S{uj)) and cry = (Ty(5(cc;) with a a 
constant value of order 1/3. 




o-j(O) 
1.0+ 



Figure 3: Chiral vortical (left) and magnetic (right) conductivities as function of the 
frequency at r = 0.008. Red doted points represent real part and thick blue line the 
imaginary conductivity. The real part of at w = is outside the range of the plot. 



4. Conclusion and Outlook 

We have derived Kubo formulas for the anomalous magnetic and vortical conductivi- 
ties ^B, as they are defined through the constitutive relations in the Landau frame 
for an anomalous current. These Kubo formulas can be used in order to perform a 
first principle field theoretical computation of the chiral vortical effect in a weakly 
coupled anomalous theory in an analogous way as it has been done for the chiral 
magnetic effect in [6]. Such result would allow us to confirm that the chiral vortical 
effect is coupling independent, as it is expected since it is anomaly induced. 

The peculiar form of the Kubo formula relies on the fact that these coefficients 
are defined for the local rest frame, so they measure the current generated in the 
presence of external background fields after subtraction of the energy flow contri- 
bution. This flow contribution is fixed once we have chosen a frame (e.g. Landau 
frame) for the definition of the fluid velocity. For this reason we have also defined 
elementary conductivities, as and ay, that measure the complete response to the 
external magnetic and gravitomagnetic field. The analytic results found for the zero 
frequency conductivities were given in (3.41) - (3.44). They reproduce the known 
holographic results for the ^ conductivities. The chiral magnetic conductivity aB 
also reproduces the weak coupling results for an axial anomalous theory, whereas the 
chiral vortical conductivity ay is a genuine new result. The frequency dependence 
of the aB shows that at very small temperatures, the conductivity drops to 1/3 of 
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its zero frequency value as soon as one moves to w 7^ 0. On the other hand, even for 
large temperature, oy decays very fast as a function of the frequency. It would be 
interesting to have a weak coupling computation of the chiral vortical conductivity in 
order to see if the strong interactions have some effect [45]. In principle, both effects 
are experimentally accessible since they induce a macroscopical effect, they modify 
the hydrodynamic description of the system. In fact, they have been argued for 
being relevant in heavy ion collisions to explain the charge asymmetry fluctuations 
observed at RHIC for non-central collisions. 

A. Second Order Action 

The general contributions to the second order action in perturbations of the gauge 
field and the metric coming from the bulk action, the Gibbons-Hawking term and 
the gravity and gauge counterterms, respectively, are 

+ \ [g^^'F^'' + 2 {g^^F^^ - g^^F^^)) h^sac] , (A.l) 

4 J r— >oo 

= I d'xV^{l {n^-g^^g^^ - 2nV^^?^^) hAsS^chnE (A.2) 

where is the normal vector to the boundary. 



+ ^ V, VV"/^^'^ - V.h^^'V^K " - \v^h^^^V,K °) ] , (A.3) 

+ F.^F^ph^^'^h'P - F, -F^M^^ + 2 {g'^'^g''^ - g^^g'"') V^a.V^ap 
+ 2hF^^V^a, - Ah^'^F^ "V^a,} . (A.4) 
From these terms we extract the A and B matrices in the boundary term (3.16). 
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B. Analytic solutions in the hydrodynamic regime 



The perturbative solutions of the system (3.24) and (3.25) up to first order in mo- 
mentum are 



hliu) = hlfiu) - 



ikKeij{u — l)a 
2(l + 4a)3/2 



2 (l + M - au'^) ArcCoth 



3 i^Vl + 4an(2an - 1 
2 + u 



B.J+ (l + 4a)3/V/i^' 



BAu) = Bi + h\u - i 



2(1 + 4a)3/2 



Bj (^6ay/T+ iau + (-2 + a(-2 + 3M))Log 



Uu (1 + 4a)'/' + 

2 + u- y/1 + 4an 



(B.l) 



2 + U + VI + iau 



■ (B.2) 



C. Fluid velocity dependence 



We have seen that in the hydrodynamic regime the velocity of the fluid in the Landau 
frame is determined modulo a P-odd term Vm ~ 0{k) that is an arbitrary function of 
the sources. In this appendix we show the independence of the transport coefficients 
on this arbitrary function, even if the correlators are velocity dependent, and also 
that these arbitrariness disappears once we correctly impose the physical boundary 
conditions on the bulk fields. 

In order to do so, we are going to solve the system at a; = 0, first order in k and 
for arbitrary value of Vm- Again, the system reduces to: 



= h'f{u) 



= BAu) 



h'iiu) 



3au B'Au) 



u 



^'^""^ BAu) -re 



(C.l) 
(C.2) 



where h\{u) = hf\u) +phf'\u) and Bi{u) = Bf\u) +pBf^{u). After imposing 
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regularity at the horizon we find the following solutions: 



B\'>{u) = C,u 



B, + AiU, 

2i{l + a^ReijAjU 



(2 - a) (1 + 4a)6 (2 - a) (1 + Aafl% 



+ (2 - af yiil + a)ArcCoth 
hf\u) = K + A, 



Vl + 4a 



+ (2 + a(2 - 3n))ArcTanh 



(C.3) 

9a(l + a)MVl + 4a (C.4) 
1 + 2aM 



VI + 4a 



i a fi;(4(l + a^u — 27a) tijAj v? 



3ia K eijBj 



2(2-a)(l + 4a)3/25 



2(2-a)(l + 4a)6 
(^((2 + a(16 + 5a))M - 6a(l + a)^^ - (2 - a)^)^! + 4aM 

-1 + 2aM 
\/l + 4a 



(C.5) 
(C.6) 



+2(2 - a)2 I^ArcCoth Vl + 4a + /(M)ArcTanh 



As we know, this is not enough to solve the boundary value problem since both of the 
two independent solutions for the metric fluctuations satisfy the regularity condition. 
However, we can use the constitutive relations to try to fix the arbitrariness. In the 
hydrodynamic description, the stress-energy tensor is given by 



T*^ = (e + Py - Ph^ 



(C.7) 



where the velocity is order p. Using the holographic dictionary, we can identify the 
coefficient of the non-normalizable mode of the asymptotic behavior of a bulk field 
with the source of the dual operator and the coefficient of the normalizable one with 
its expectation value. Therefore, we can write the metric fiuctuation close to the 
boundary as 



h\{u)^h\ + Tlu\ 



(C. 



so using the hydrodynamic result, we can do the identification order by order in 
momentum, in such a way that the velocity piece of the energy tensor fixes the 
horizon value of h^l'\ Doing so, the asymptotic behavior of each order becomes 



ph. 



hf\u) ^ K{1 - Pu' 
(1), 



u) ~ 



-{e + Pyu" 



(C.9) 



We can proceed to construct the matrix of correlators for arbitrary value of the veloc- 
ity as explained in section 3. Now, all the correlators pick contributions proportional 
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to the velocity. In a compact way, the retarded propagators read 



~ ttGL [ 8(1 + a) 2L2 QB^ ' ' ^ ' 



_ 'H I 3iakK y/Sar-n v^u-'h i^l'j i /'r'li'i 
G.,,-^'-^^, (C.12) 





\^4(l + a) 


4 (1 


+ a) dvi 






4 ^ 


+ 




where i,j = x,z. It is straightforward to prove that applying definitions (1.3) and 
(1.2) for the chiral vortical and magnetic conductivities, the result is independent of 
the velocity and coincides with (3.44) and (3.43) as expected. Setting the velocities 
to zero, the correlators coincide with those presented in [43]. 

If we now impose the correct zero frequency 'infalling' condition to the fields h], 
i. e. vanishing at the horizon, the velocities are not arbitrary anymore, but are given 
in terms of the boundary sources, 

iaReij{2hj + 3Bj)k 
16(1 + a) 

Of course, substituting them in the Green functions given above, the antisymmetric 
correlation matrix spanned by (3.34) - (3.38) is recovered. 
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